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A model for A = 3 antinuclei production in proton-nucleus 
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Abstract. A simple coalescence model based on the same diagrammatic approach of antimatter production 
in hadronic collisions as used previously for antideuterons is used here for the hadroproduction of mass 3 
antinuclei. It is shown that the model is able to reproduce the existing experimental data on t and '^He 
production without any additional parameter. 

PACS. 24.10.-i Nuclear-reaction models and methods 



1 Introduction 

The increasing interest in the study of production of light 
antinuclei in proton-proton and proton-nucleus collisions 
is motivated by the presence of anti-nuclei in cosmic rays 
which has potentially important implications on the matter- 
antimatter asymmetry of the universe. From this point of 
view, it is important to determine the amount of anti- 
matter which can be produced in the galaxy through the 
interaction of high-energy protons with the interstellar 
gas. A new generation of experiments (AMS PAMELA 
0) should be able to measure the flux of anti-matter in a 

near future. 

] The calculations of the t and -^He production cross sec- 
. tions reported here are based on the same diagrammatic 
' approach to the coalescence model as used recently |31 to 
describe the d production in proton-proton and proton- 
nucleus collisions. 

The coalescence model 0] is based on the simple hy- 
pothesis that the nucleons, produced during the collision 
of a beam and a target, fuse into light nuclei whenever 
the momentum of their relative motion is smaller than a 
coalescence radius po in the momentum space, which is 
a free parameter of the model, usually fit to the experi- 
mental data (see ^ for example). A simple diagrammatic 
approach to the coalescence model developed in [H| pro- 
vided a microscopic basis to the model. In this approach, 
the parameter po is expressed in terms of the slope pa- 
rameter of the inclusive nucleon production spectrum and 
of the wave function of the produced nucleus. 

This diagrammatic approach has been generalized in 
|3] to antideuteron production by taking into account thresh- 
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old effects and the anisotropy of the angular distributions. 
This approach can reproduce most existing data with- 
out any additional parameter in energy domains where 
the inclusive antiproton production cross sections are well 
known. 

This article reports on the application of this approach 
to the production of ^ 3 antinuclei. It is the first micro- 
scopic calculation of this cross section to the knowledge 
of the authors. In [7], the '^He production cross section in 
proton-proton collisions was calculated using the standard 
coalescence model, with the parameter po taken from the 
d production data. 

Unfortunately, the experimental data required to be 
compare to the calculations are limited. Only two sets 
of experiments have measured the production of mass 
3 antinuclei in proton-nucleus collisions, t and '^He were 
discovered at IHEP (Serpukhov), with one experimental 
points measured for t and one for ■^He |HE1, while in the 
CERN experiment (SPS, WA 33) [TIlTTl . four experimen- 
tal points were measured for t and eight for ■^He. For these 
latter data however, the t and -^He production cross sec- 
tions were measured with respect to the pion production 
cross section at the same momentum. This requires the 
corresponding experimental values of the pion production 
cros s sec tion to be known to extract the values of the t 
and ^He production cross sections. 

The article is organized as follows. The main ideas of 
the theoretical approach are described in section 2. The 
formalism is generalized to the case of ^ = 3 antinuclei 
production in section 3. Section 4 is devoted to the results 
and the comparison to the experimental data. A brief sum- 
mary of the work is provided before the work is concluded 
in the last section. 
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2 Diagrammatic approach to the coalescence 
model 

The main ideas of the diagrammatic approach of the co- 
alescence model for nuclear fragment production arc re- 
minded here for the reader's convenience [H]. The sim- 
plest Feynman diagram of Fig. Q corresponding to fusion 
of three nucleons is considered as a basis for the coales- 
cence model. Here the symbol / designates the state of all 
particles but nucleons 1, 2 and 3 which form the tritium or 
the helium 3 nucleus produced in the final state (specified 
by the t symbol on the graph). 
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Fig. 1. The simplest Feynman diagram corresponding to coa- 
lescence of three nucleons into a tritium or ^He. 



The physical picture behind this diagram is quite sim- 
ple: the nucleons produced in a collision (block A) are 
slightly virtual and can fuse without any further inter- 
action with the nuclear field. This diagram is not the 
only possible contribution to the full transition amplitude. 
However, mutual cancellations of a number of other con- 
tributing diagrams result in the diagram of Fig. ^ being 
dominant This diagram can be calculated using the 
technique developed in |13j . 

The probability for three nucleon coalescence is given 

by: 



since they cancel in further calculations. Using the con- 
ventional graph technique |13j . the expression for M can 
be written in the form: 
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where M(i 2,3-»t) is the vertex of coalescence of 1,2,3 into 
t (proportional to the three-nucleon wave function in the 
momentum space in the nonrelativistic approximation). 
TOp the nucleon mass, the three fractions being the in- 
dividual nucleon propagators of 1, 2 and 3. The inte- 
grals have to be performed over energies and momenta of 
the (virtual) particles. The delta functions ensure energy- 
momentum conservation at the t vertex, pt = {'PnEt)^ 
with Pt = Pi -|- P2 + P3 being the momentum of t, Et 
its energy. The dependence of the amplitude Ma on its 
variables (the particle momenta) is also needed explicitly 
for the calculations. In lack of a reliable theoretical form, 
this can be done in a "minimal" way, by using empirical 
shapes. The inclusive nucleon spectra usually have a de- 
creasing form which can be approximated by a Gaussian 
function in the center of mass frame: 



^^P$^«exp(-p^/g^) 



(4) 



where Q defines the slope parameter of the momentum 
distribution. Accordingly, the amplitude Ma can be writ- 
ten in the following way: 
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where rrit in the mass of the t fragment and Et its energy in 
the (beam)nucleon-(target)nucleon center of mass system 
of the colliding nuclei. The probability for three nucleon 
production is then: 
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where Ma is the amplitude corresponding to the block A, 
i.e., accounting for the inclusive production of nucleons 1, 
2 and 3 and other particles in the final state /. To avoid 
cumbersome expressions in equations and Q , the fac- 
tors corresponding to colliding nuclei have been omitted 



Pt = Pi +P2 + P3, 
P = — (Pl - P2) , 
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Assuming a statistical independence in the three nucleon 
production process, the inclusive production cross section 
can be written as the product of the three independent 
probabilities: 
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where Umei is the total reaction cross-section of the col- 
liding particles. 
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After integration of (j^J, taking into account (jSJ and 
l|7|l- and dividing by the incident particle flux, the t pro- 
duction cross section takes the form: 
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with Pi = P2 = P3, Pt = 3pi and 
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Where tf't (p,q) fx Afi23^t is the wave function of the t 
(or "^He) system normalized by the condition 
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(10) 



The factor 1/2, accounts for nucleons and ^ = 3 nuclei 
spins, is included in ((HJ. The three- nuclcon wave function 
is needed at sufficiently large momenta to compute the 
amplitude. The wave function of has been used (see 
appendix for discussion). 

The structure of (jSj) is the same as that of the coales- 
cence model and the S integral in |H1 can be straightfor- 
wardly related to the coalescence momentum: 
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Thus, in the approach based on the diagram of Fig. ^ 
and within the approximations made above, the coales- 
cence momentum po is not an adjustable parameter any- 
more, but it is determined by the inclusive proton spec- 
trum and by the trinuclcon wave function. Note that in 
that case, po depends on the momentum distribution and 
should thus be energy and system dependent. 



3 Application to three-antinuclei production 

In order to generalize the diagrammatic approach of the 
coalescence model to the production of A = 3 antinu- 
clei (noted t further below), two effects have to be taken 
into account: the anisotropy of angular distributions and 
the threshold effects Isotropic angular dependence are 
frequently assumed in nonrelativistic collisions. However, 
in relativistic collisions, the momentum distributions are 
strongly anisotropic and the low energy approximation 
cannot be used. To take this into account, formula ((SJ 
can be easily generalized to any angular dependence. As- 
suming the inclusive nucleon production cross section to 
be given by the (parameterized) amplitude AIi (pi): 



i^i^ = iA/i(Pi)r 

dpi 



(12) 




Fig. 2. Dependence of the threshold factor R{x), with x as 
defined in the text. 



The cross section for t production can then be written (see 
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This model could be practically used directly to describe 
the production of t. The production threshold of the an- 
tiparticlc has to be taken into account however in the cross 
section calculation. The same procedure to evaluate the 
cross section near threshold as that used in [2j will be 
applied here. In nucleon-nucleon collisions, the main re- 
action producing a 1 particle is NN t + 5A^. Near the 
threshold of this reaction, the energy dependence of the 
t production cross section is mostly governed by the five 

nucleons phase space: ^ (^\/ s + — 2^/sEi; 5mp^ , 



d^aj 
dpi 



s + nil - 2\/s£'t; 5r 



(14) 



The phase space <!> for n particles with masses, momenta 
and energies, m^, Pj, and Ei respectively, is defined in the 
usual way (in the center of mass) 



1 d^p, ,3 



n 



\ (27r)3 2E, 



It was calculated here by using the standard CERN library 
program (W515, subroutine GENBOD) |T3]. is the 
total energy of the n particles in the center of mass system. 

A phenomenological correction factor R can thus be 
introduced in formulae 113|l which then reads: 



R{x) 



<P (cc; 5mp) 
<Z'(x;5 X 0)' 



(15) 
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where x = \J s + — and where the denominator 

contains the high energy hmit of the phase space to ensure 
R to be dimensionless and to do not change the value of 
the cross section out of the space phase boundary. The 
hmits of R are thus: 



i? ^ 0, E,^E^' 

R ^ 1, y/s CX) 



s + ml — {brripY 



mt 



the expression \/ s + ' 
-Et 



Ifp|«(Vi-i5,)2, 
can be replaced by ^/s — E-. This same approximation was 
made in The functional dependence of R{x) is shown 
in Fig. El 



4 Results on antinuclei production data 

4.1 Status of the data 

This section is introduced with a brief overview of the 
current experimental situation on the antinuclei produc- 
tion relevant to the present study, i.e., in proton-proton 
and proton-nucleus collisions. The antinuclei production 
in ion-ion collisions will be quoted only for completeness. 

— As mentioned in the introduction, the experimental 
data on the production of mass 3 antinuclei are ex- 
tremely scarce, and much less informative than that on 
antideuteron production, with only two experiments or 
sets of experiments reporting on mass 3 antinuclei pro- 
duction in proton- nucleus collisions [SlIUlllOlllH . Note 
that there are no experimental data available on the 
production of these antinuclei in proton-proton colli- 
sions. The production of '^He has been observed re- 
cently in various heavy ion studies like Pb + Pb colli- 
sions at ultra relativistic incident energies )16| . These 
data are out of the scope of the present work. They 
will not be discussed here (see |3| for a discussion). 

— Coalescence calculations require the antiproton pro- 
duction cross section to be known for antiproton mo- 
menta equal to approximately one third of the A = i 
antinuclei momenta. Unfortunately, in most experi- 
ments the differential cross sections for antiproton and 
^ = 3 antinuclei productions were not measured at 
this momentum. The p cross section thus had to be ex- 
trapolated to the appropriate kinematical region when 
no other data were available, which of course, intro- 
duces additional uncertainty in the calculations. 

The three nucleon wave functions needed in the calcu- 
lations are much less well know than the deuteron wave 
function. In addition, the same wave function will be used 
for ^He and t nuclei (see appendix). The inaccuracy on 
the tri-nucleon wave functions is thus another source of 
uncertainty. 

The total reaction cross-section used in the calcula- 
tions was described by means of the parameterization pro- 
posed in |17| . 



4.2 Proton Aluminium collision data at 70GeV/c 

The antinuclei produced in the Serpukhov experiments [HI 
El were obtained from a 70 GeV/c proton beam incident 
on an aluminium target at 27 mrad scattering angles and 
20 GeV/c for ^c, and and 25 GeV/c for t. The inclusive 
antiproton cross sections were available from and [HI] 
in the same kinematical conditions. 

In Fig. El the p cross section data from |H1 are com- 
pared with the results of fits using a functional form |20| . 
The solid curve corresponds to a fit of a large sample of 
p + A—tp data from 12 up to 400 GeV incident energies 
not including those from reference which were found 
not to be compatible with the other sets of data (21] . The 
calculated values are in fair agreement with the two lowest 
momentum data points (which were obtained by extrap- 
olation from measurements at other angles). They over- 
estimate the other data points by a factor of 2 to 4. The 
dotted curve is a renormalization of the solid curve by 
a factor sa 2.5 to fit these latter points, while the dashed 
curve corresponds to the fit of the single set of data points 
shown on the figure, which parameters however give quite 
poor agreement with the other sets of data |20| . 
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o 




50 60 

Piab (GeV/c) 

Fig. 3. Inclusive differential cross section for antiproton pro- 
duction in p-|-Al collisions as a function of the total momentum 
in the laboratory frame from [SI, compared to calculated values 
as discussed in the text. 



It must be emphasized that the low momentum region, 
say piab < 10 GeV/c, which is the useful region for the 
coalescence calculations, with Pp « p-/3 is particularly 
important here, with unfortunately no data point from 
direct measurement available over the relevant range. 

Fig. El compares the calculations for the A = 3 produc- 
tion cross section for the three parametcrizations shown 
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P..(GeV/c) 

Fig. 4. The inclusive diflterential cross-section for t |^ (open 
circles) and ^He [5| (black circles) production by 70 GeV pro- 
tons on Al target as a function of the total momentum in the 
laboratory frame compared to calculations using the micro- 
scopic model of coalescence and the same three different pa- 
rameterizations of the antiproton production cross-section as 
shown on figure |21 with the same graphical conventions. 



parameter Ci C2 C3 C4 
value 0.94 1.88 7.05 1.69 



Table 1. Values of the parameters of relation^] obtained by 
fitting the n~ production cross sections for 200 and 300 GeV/c 
protons on Beryllium. 



experiment. The measured distributions have been fit by 
means of the following funetional form, inspired from ref 

m 

(tt-) ^ Cia,„ (1 - xf' e-^^-e-^=f- , (16) 

where x = E~ /E^^^ {E* is the total energy of the inclu- 
sive partiele in the center of mass frame, <Jin is the total 
reaction cross section for the system in collision, ^/s is the 
total energy of the system and the transverse momen- 
tum of the emitted particle. The values of the parameters 
obtained are given in TableQland the results of this fit are 
presented in Fig. The parameterization H16() has been 



on Fig. with the experimental data. The calculations 
using the global fit renormalized to the 70 GeV/c data 
(see Fig. give by far the best agreement with the t data 
(dotted curve). The other two sets of p cross section pa- 
rameters overestimate the data by a sound order of magni- 
tude. This is apparently consistent with the larger p cross 
section predicted by these two sets of parameters for low 
p momentum region to which the t cross section is most 
sensitive. The factor of about 2 between the p cross sec- 
tions predicted by the two groups of parameters translates 
into a factor of about 10 for the t cross section because 
of the approximately cubic dependence of the latter on 
the p cross section. However it is somewhat puzzling that 
this agreement is obtained with parameters which are not 
consistent with the whole body of p data [20] . 



4.3 Proton beryllium collision data at 200 GeV/c 

In the CERN experiments |10llllj . p, t and •^He were pro- 
duced in proton-beryllium collisions at 200, 210, and 240 
GeV/c and detected at degree scattering angle [TT] . 
while p were measured at 200 GeV/c ^01 on the same 
targets. For these data however, the production cross sec- 
tions were measured as the ratios to the tt" production 
cross sections at the same momentum. The knowledge of 
the corresponding experimental 7r~ production cross sec- 
tion, or a good parameterisation of the latter, is thus re- 
quired in order to allow the values of the p, t and -^He 
production cross sections to be calculated. 

Fortunately, the p + Be — > tt" -I- X cross section has 
been measured at 200 and 300 GeV/c incident momentum 
in |2J in similar kinematical conditions as in the CERN 



p-bB-200-JOO GeV/c 




40 60 80 100 120 140 '60 180 

p„ (GeV/c) 



Fig. 5. Experimental inclusive differential cross-section for 7r~ 
production in p -|- Be collisions |21) (symbols) at 200 GeV/c 
(full circles) and 300 GeV/c (full triangles) compared with the 
functional form 116! 



used to extract the experimental p production cross sec- 
tions IT^ . The resulting cross section values are compared 
in Fig. El with the results of the fit of a functional form to 
a large sample oip + A^p-\-X data from 12 GeV/c up 
to 400 GeV/c incident momenta It is seen that the 
data points derived previously and the calculated values 
are in fair agreement. This consistency gives confidence 
to the following steps of the analysis for the evaluation 
of the t and ^He production cross sections. In Fig. [T] 
the t and ■^He production cross sections are compared to 
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p-AI-Be-20Q GeV/c 



40 45 50 
P» (GeV/c) 



Fig. 6. Experimental inclusive differential cross-sections for p 
production in p + Al collision (circles), and in p + Be collision 
(xlO~^, triangles), [18] compared with the results of fits using 
a functional form (curves) |20|. 



calculations using the microscopic model of coalescence. 
The agreement between experimental and calculated val- 
ues varies from poor to good. On the average however data 
and calculations are within one order of magnitude. This 
result should be considered as a success in account of the 
numerous sources of uncertainties of the calculations and 
of the limited accuracy of the measurements. Note also 
that refs ^Uj and ^Jj report experimental values in dis- 
agreement by a factor of 2. Furthermore the t and ^He 
production cross section, measured at the same momen- 
tum, should be in principly close to each other (this fact is 
clearly seen in the same experiment for t and ^He produc- 
tion) whereas, in this experiment, they are quite different. 





p+Al-200 GeV/c 




• 



Pi_^ (QeV/c) 



p+Be-200 GeV/c 




■ 



10 15 20 25 30 35 40 45 50 

p,_JGeV/c) 




5 Conclusion 

It has been shown in this work that the diagrammatic 
approach to the coalescence model developed previously 
previously can successfully account for the mass 3 antin- 
uclei production cross section in proton-nucleus collisions 
over wide kinematical conditions without any additional 
parameter. These calculations require a good knowledge of 
the antiproton production cross-section and of the three- 
nucleon wave function. Th ese results would be further 
used to calculate t and '^He flux in cosmic rays. 



Appendix 

In this appendix, we briefly remind how the wave func- 
tions of the trinucleon is written in |14| . while in this pa- 
per slightly different definitions have been used. A useful 
analytical parameterization of the bound trinucleon wave 
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Fig. 7. Experimental inclusive differential cross-section for 1 
(full circles) and ^He (open circles) production in p+Al and p+ 
Be collisions, compared to calculations using the microscopic 
model of coalescence (solid line) . 



function is obtained from solving the Faddeev equation 
with the Reid soft-core potential. 
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The total wave function of the triton is written as a 
sum of 3 Faddeev components 

3 
i=l 

In (|5|l- 113|l . we make use of only one Faddeev component 
^t, while due the exchange symmetry of two-nucleons in 
the triton, all the three Faddeev components are identical. 
Furthermore, Faddeev components are decomposed in 
terms of their partial wave components with respect to 
the spin-isospin and angular momentum basis '?!>c((p,q) ■ 

with, if Pi (« = 1, 2, 3) is the nucleon momenta 

P = ^ (Pi - P2) , q = ^ (Pi + P2 - 2p3) . 
The following normalization is used 



l>ft(p,q) 



■ (i^pd^q 



Note that, in these expressions, the definition of p and 
the normalization differ from ((HJ and H10() . tf^t is a sum of 
the partial wave state a which is a label for the following 
physical quantities: 

— L, the angular momentum of the pair of nucleons (1-2). 

— I, the angular momentum of nucleon 3 according to the 
center of the mass of the pair of nucleons (1-2). 

— £, the total angular momentum of the triton. 

— s, the spin of the pair of nucleons (1-2). 

— iS, the total spin of the triton. 

~ T, the isospin of the pair of nucleons (1-2). 

Only two components label a were taken into account, 
a = 1,2. 



— For a 

— For a 



1, L : 

2, L 



C 
C 



0. 



0, s = 1, 5 = 1/2 and T 
s = 0, 5 = 1/2 and T = 1. 

Of course, the fact to consider only two partial wave state 
is an approximation which gives the probability of 89.25% 
of trinucleon being in the partial wave state a. In |14j . the 
parameterization for ipa {p, q) is given by 
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^„ (p, q) = p^p' + nl^) ' n + ^IraY 



6 6 

EE 



-[U{p^+i4){i^+^^y 

with i7pi, ilqm, Hi, Vj and Cij all depending on the partial 
wave label a. The numerical values of these coefficients can 
be found in |14j . 

^He{ppn) and t{pnn) are considered to have the same 
wave function. Altough, because of the presence of the 
Coulomb interaction, these two wave functions are slightly 
different, this difference is negligible compared to the other 
uncertainties of present calculations. 



